屁 虫 学 报 http : Nwww. insect. org. cn 
ACTA ENTOMOLOGICA SINICA doi: 10.16380/j. kexb. 2019. 05. 003 











花椒 窄 吉本 转录 组 及 化 学 感受 相关 基因 的 分 析 


杨 平 ， 谢 寿 安 ”， 项 雪 芳 , 车 显 荣 , 王 延 来 , 吕 淑 杰 


(西北 农林 科技 大 学 林学 院 , 陕西 杨陵 712100) 











摘要 :【 目 的】 建立 花椒 窒 吉 丁 Agrilus zanthoxylumi 转录 组 数据 库 , 挖 据 其 基因 组 数据 。 【方法 】 采 
用 高 通 量 测序 平台 Ilumina NovaSeq6000 对 花椒 窄 吉 丁 成 虫 不 同 组 织 ( 触 角 、 头 、 胸 、 腹 、 足 、 起 ) 进 
行 转录 组 测序 序列 组 装 ,使 用 BLAST 软件 将 花椒 窄 吉 丁 unigenes 序列 与 公共 数据 库 进行 比 对 ， 
BLAST 同 源 性 搜索 的 方法 从 中 筛选 出 与 其 化 学 感受 相关 的 基因 ,并 与 其 他 已 研究 发 表 的 鞘翅 目 昆 
来 化 学 感受 相关 基因 的 核酸 序列 比 对 ,进行 系统 发 育 分 析 。 利 用 RPKM 值 对 这 些 基因 在 花椒 窄 吉 
丁 成 蛙 不 同 组 织 中 的 表达 量 进行 分 析 。【 结果 】 经 测序 及 序列 拼接 后 共 得 到 80 320 条 unigenes 和 
169 398 条 contigs ,其 G +C 比例 分 别 是 37.63% 和 39.18% ,平均 长 度 分 别 为 828.75 和 1 084. 33 
bp; unigenes 长 度 主要 分 布 在 200 ~400 bp 的 共有 37 374 条 , 占 全 部 unigenes 的 46.53%。 在 NR 
数据 库 中 注释 的 与 花椒 罕 吉 丁 转录 组 相似 基因 所 属 物种 分 布 为 赤 拟 谷 资 Tribolium castaneum 所 占 
比例 26. 80% ,其 后 依次 是 小 家 鼠 Mus musculus(13.87% ) , 驻 犀 金龟 Oryctes borbonicus(5.26% ) , 山 
松 甲 于 Dendroctonus ponderosae(4.10% ) , 黑 Lasius niger(2.36% ) ,其 他 物种 所 占 比 例 是 47.61% 。 
在 GO 数据 库 中 比 对 到 27 488 个 unigenes , 可 分 为 分 子 功能 、 细 胞 组 分 和 生物 学 进程 三 大 类 共 59 
个 过 程 ;利用 KOG 数据 库 功 能 注释 分 为 25 类 ,注释 到 unigenes 最 多 功能 类 别 的 是 普通 功能 , 共 
4 666 个 ,最 少 的 是 与 细胞 移动 相关 , 共 40 个 ;将 80 320 个 unigenes 映射 到 KEGG 数据 库 中 ,21 104 
条 unigens 注释 到 代谢 通路 35 个 , 占 26.27% ,注释 最 多 的 代谢 途径 是 转录 途径 , 共 涉 及 到 2 037 条 
unigenes。 从 NR 数据 库 中 共 注 释 到 8 个 气味 受 体 (odorant receptor，OR ) 基因 、7 个 离子 型 受 体 
(ionotropic receptor, IR) 基 因 和 6 个 味觉 受 体 (gustatory receptor，GR ) 基 因 ,各 组 织 中 的 基因 表达 量 
分 析 发 现 ,这 3 类 化 学 感受 相关 基因 在 成 来 触角 中 的 表达 量 显 著 高 于 其 他 组 织 中 的 表达 量 , 上 只 、 雄 
触角 中 表达 量 差异 最 显著 的 为 AzanORI ,而 在 肉 、 雄 成 虫 足 中 这 3 类 化 学 感受 相关 基因 表达 量 没有 
差异 。【 结论 】 本 研究 首次 获得 了 花椒 窄 吉 丁 转录 组 数据 ,为 进一步 研究 花椒 窄 吉 丁 的 基因 功能 葛 
定 了 分 子 基础 。 
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Analysis of the transcriptome and chemoreception-related genes of Agrilus 


zanthoxylumi ( Coleoptera: Buprestidae ) 
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( College of Forestry, Northwest A&F University, Yangling, Shaanxi 712100, China) 

Abstract: [Aim) To establish the transcriptome database of the Agrilus zanthoxylumi, and to obtain its 
genomic data. 【 Methods) The transcriptome of the different tissues (antennae, head, thorax, abdomen, 
leg, wing) of A. zanthoxylumi adults was sequenced using an Illumina NovaSeq 6000 platform, and the 


unigene sequences were compared with the homologous sequences in public databases by BLAST. 
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Chemoreception-related genes were Screened out by BLAST homology search, and their nucleotide 
sequences were compared with those of other reported chemoreception-related genes of Coleoptera insects 
for phylogenetic analysis. The expression levels of these genes in different tissues of A. zanthoxylumi 
adults were analyzed based on the RPKM value. 【Results】 After sequencing and sequence splicing, a 
total of 80 320 genes and 169 398 contigs were obtained. The G + C proportions of contigs and unigenes 
were 37.63% and 39. 18% ，respectively ，and the average length of contigs and the average length of 
unigenes were 1 084.33 and 828. 73 bp, respectively. There were 37 374 unigenes with a length 
distribution mainly in 200 - 400 bp, accounting for 46. 53% of the total unigenes. Origin species 
distribution of similar genes with unigenes of A. zanthoxylumi annotated in the NR database was Tribolium 
castaneum (26. 80%), Mus musculus (13. 87%), Oryctes borbonicus (5. 26% ), Dendroctonus 
ponderosae (4. 10% ), Lasius niger (2. 36% ) and other species (47. 61% ) There were 27 488 
unigenes aligned to the GO database, which can be divided into three functional categories including 
molecular function cellular component and biological process in a total of 59 processes. According to the 
gene annotation against the KOG database, these unigenes can be divided into 25 categories. The most 
annotated functional category was the common function (4 666 unigenes), while the least annotated was 
cell movement-related (only 40 unigenes). When 80 320 unigenes were mapped to the KEGG database, 
21 104 unigenes were grouped into 35 metabolic pathways, accounting for 26. 27% of the total. The most 
annotated metabolic pathway was the transcription pathway (2 037 unigenes). By gene annotation against 
the NR database, we found eight odorant receptor (OR) genes, seven ionotropic receptor (IR) genes 
and six gustatory receptor ( CR ) genes. Tissue expression profiles showed that the three kinds of 
chemoreception-related genes (OR, IR and GR genes ) had significantly higher expression levels in the 
antennae than in other tissues of adults and the expression level of AzanOR!1 in the antennae showed the 
most significant difference between female and male, while the expression levels of the three kinds of 
chemoreception-related genes in the legs of adults showed no difference between female and male. 
【 Conclusion】This study obtained the transcriptome data of A. zanthoxylumi. The results provide a 
molecular foundation for further studying the gene function in A. zanthoxylumi. 
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花椒 窗 吉 丁 Agrilus zanthoxylumi 隶属 峭 怒 日 
( Coleoptera ) 吉 丁 虫 科 (Buprestidae ) 罕 吉 丁 属 
hgrilus( 孙 益 知 和 李 忠 锋 ，1988 ) ,一 化 性 昆虫 , 均 为 
一 年 发 生 1 代 , 主 要 分 布 在 我 国 北方 花椒 产 区 ,比如 
陕西 甘肃、 山西 ` 山 东 等 。 目 前 研究 发 现 花椒 窄 吉 
丁 仅 危害 花椒 ,为 单 食性 昆虫 ( 吴 海 , 2006) 。 幼 虫 
在 木质 部 或 树 皮 隧道 中 越冬 ,成 虫 以 叶片 为 食 , 对 花 
概 叶 造成 缺 刻 和 空洞 (和 崔 艳 鱼 ，1990 ) ,受害 花椒 树 
皮 的 主要 症状 是 会 大 量 流 胶 、 软 化 腐烂 等 , 树 皮 干 
枯 龟 裂 ,严重 会 导致 伦 椒 树 体 死亡 ( 何 凯 ，1991 ) 。 
现 阶 段 对 花椒 罕 吉 丁 的 研究 主要 集中 于 化 学 防治 和 
生物 学 特性 两 方面 ( 康 克 功 等 ,2007 ) ,利用 天 敌对 
其 进行 防 控 也 有 少量 报道 ( 张 伟 ,2009 ) ,而 长 期 使 
用 化 学 药剂 防治 的 策略 会 使 害虫 产生 3R 问题 ( 张 
美 滩 等 , 2014) ,会 导致 运用 化 学 防治 也 很 难 奏效 ， 
还 会 影响 花椒 的 产量 及 其 商品 价值 ,给 花椒 产业 造 



































成 了 重大 的 经 济 损失 。 生 物 防治 作为 应 用 较为 广泛 
的 无 公害 防治 技术 之 一 ( 陈 学 新 等 , 2013) 。 随 着 我 
国 绿色 无 公害 生物 防治 技术 的 发 展 , 探 究 农 林业 有 
害 昆 虫 的 高 效 防 控 ,基因 组 信息 不 全 已 成 为 吸 待 解 
决 的 重要 问题 之 一 。 然 而 ,关于 花椒 罕 吉 丁 转录 组 
言 息 方 面 的 研究 还 是 空白 。 

选用 高 通 量 转录 组 测序 技术 测定 花椒 窄 吉 丁 转 
录 组 信息 的 原因 ,主要 是 高 通 量 转录 组 测序 技术 拥 
有 精确 到 单个 核 苷 酸 的 分 辨 率 , 且 能 对 样品 任意 时 
间 点 或 任意 条 件 下 的 转录 组 进行 测序 ,还 能 够 动态 
反映 基因 转录 水 平 ( 周 华 等 , 2012) ,同时 鉴定 和 定 
量 稀 有 转录 本 和 正常 转录 本 ,提供 样品 特异 的 转录 
本 序列 结构 信息 和 确定 基因 序列 和 转录 本 ,在 没有 
参考 基因 的 条 件 下 也 可 以 对 某 一 物种 的 基因 组 进行 
细致 全 面 的 分 析 ,为 基因 组 图 谱 尚 未 完成 或 数据 信 
上 息 匮 乏 的 物种 的 基因 组 和 转录 组 研究 提供 技术 支持 
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等 优势 而 在 昆虫 转录 组 学 研究 上 得 到 广泛 使 用 ( 贺 
华 良 等 , 2012)。 从 2008 年 至 今 ,人 研究 人 员 已 采用 
新 一 代 测 序 技术 对 7 个 目的 68 种 昆虫 进行 了 转录 
组 测序 ,其 中 由 我 国 研 究 人 员 完 成 6 个 目的 22 种 昆 
虫 的 转录 组 测序 (杨帆 等 , 2014) 。 目 前 ,昆虫 转录 
组 学 研究 生物 体 的 细胞 或 组 织 在 特定 状态 下 基因 组 
所 转录 的 全 部 mRNA , 其 反映 了 基因 在 不 同 生命 阶 
段 . 生 理 状态 .组织 类 型 以 及 环境 条 件 下 的 表达 情况 
(Yang et al., 2014)。 对 这 些 基因 进行 深入 挖 气 、 基 
因 表达 分 析 后 ,可 为 全 面 揭示 昆虫 生命 活动 中 相关 
基因 功能 .进化 以 及 昆虫 与 其 他 生物 相互 作用 等 黄 
定 基 础 。 

本 研究 采用 Ilumina NovaSeq6000 测序 平台 对 
花椒 窗 吉 丁 成 虫 转录 组 测序 、 拼 接 组 装 得 到 
unigene。 通 过 unigene 与 数据 库 比 对 ,获得 基因 功 
能 注释 、 分 类 和 代谢 途径 分 析 , 旨 在 获得 更 多 畏 翅 目 
非 模式 生物 花椒 罕 吉 丁 的 转录 本 和 更 为 全 面 的 转录 
组 信息 ,为 生物 控制 花椒 罕 吉 丁 提供 重要 的 理论 依 
据 和 相关 参考 。 于 新 一 代 测序 技术 的 昆虫 转录 组 学 
研究 ,为 昆虫 分 子 标记 的 开发 提供 了 丰富 的 遗传 学 
资源 ( Gayral et al., 2011 ) 。 袁 丽 芳 (2016 ) 利 用 固 相 
微 蔡 取 和 GC-MS 分 析 了 花椒 健康 植株 和 受 花 椒 罕 
吉 丁 危害 植株 不 同 部 位 的 主要 挥发 性 化 合 物 ,采用 
触角 电位 (EAG) 和 了 型 嗅觉 仪 测定 了 花椒 窗 吉 
对 寄主 植物 挥发 物 的 电 生 理 和 行为 反应 。 但 花椒 罕 
吉 丁 是 如 何 感 知 寄主 气味 物质 并 准确 定位 寄主 ,在 















































1.2 总 RNA 的 提取 与 质量 检测 

采用 Trizol 法 分 别提 取 花 椒 罕 吉 丁 成 虫 不 同 组 
织 (触角 \ 头 、 胸 、 腹 、 足 .这 ) 的 总 RNA ,每 个 样品 设 
置 3 个 生物 学 重复 , 1% 琼脂 糖 凝 胶 电 泳 检测 总 
RNA 的 质量 。 饥 奥 K5500 分 光 光 度 计 ( 凯 奥 , 北 
京 ) ;检测 RNA 样品 的 纯度 (0DzwoA0Dzwo 值 ) ,根据 
安捷伦 2100 RNA Nano 6000 Assay Kit ( Agilent 
Technologies, CA, 美国 ) 检 测 RNA 样品 的 完整 性 和 
浓度 。 
1.3 文库 构建 与 测序 

1.2 节 提 取 的 总 RNA 样本 检测 合格 后 ,用 带 有 
Oligo( dT) 的 磁 珠 富 集 mRNA ,向 得 到 的 mRNA 中 加 
入 片段 缓冲 液 使 其 片段 化 成 为 短片 段 ,再 以 片段 化 
后 的 mRNA 为 模板 ,用 6 碱 基 随 机 引物 合成 cDNA 
第 1 链 , 并 加 入 缓冲 液 、 4NTPs、RNaseH 和 DNA 
Polymerase | 继续 合成 cDNA 第 2 链 , 经 过 
QIAQuickPCR 试剂 盒 纯 化 并 加 EB 缓冲 液 洗 脱 。 洗 
脱 纯化 后 的 双 链 cDNA 再 进行 末端 修复 、 加 碱 基 A、 
加 测序 接头 处 理 , 然 后 经 琼脂 糖 凝 胶 电泳 回收 目的 
大 小 片段 并 进行 PCR 扩 增 ,从 而 完成 整个 文库 制备 
工作 。 构 建 好 的 文库 用 Ilumina NovaSeq6000 平台 
进行 测序 ,测序 读 长 为 PE150。 
1.4 数据 过 滤 与 拼接 预 组 装 

测序 完成 后 获得 原始 reads ,去 除 带 接 头 .ploy-N 
和 低 质 量 的 reads, 获 得 clean reads。 每 个 碱 基 测 序 
错误 率 是 由 测序 Phred 分 值 (Phred score ，Qphred ) 





























寻找 寄主 植物 过 程 中 化 学 感受 相关 基因 是 如 何 发 挥 
作用 的 ,目前 还 没有 相关 报道 。 为 此 ,本 人 研究 通过 转 
录 组 数据 中 基因 功能 注释 及 分 析 挖 掘 到 花椒 罕 吉 
的 大 量化 学 感受 相关 基因 ,为 进一步 从 分 子 水 平 开 
展 花 椒 罕 吉 丁 防治 新 策略 的 研究 莫 定 理论 基础 。 




















1 材料 与 方法 


1.1 供 试 昆虫 

花椒 罕 吉 丁 成 虫 采 集 于 陕西 省 西安 市 蓝田 县 普 
化 镇 周 董 村 花椒 园 ,采集 时 间 为 2018 年 5 -7 月 。 
背 助 捕虫 工具 寻找 在 叶片 ` 校 于 上 栖息 的 成 虫 迅 速 














通过 公式 转化 得 到 , 而 Phred 分 值 是 在 碱 基 识 别 
(base calling) 过 程 中 ,通过 一 种 预测 碱 基 判 别 发 生 
错误 概率 模型 计算 Q20，Q30 和 GC 含量 和 重复 序 
列 水 平 。 采 用 Trinity 软件 对 clean reads 进行 拼接 、 
过 滤 和 组 装 后 获得 高 质量 的 unigene。 由 于 该 样品 
昆虫 花椒 窄 吉 丁 为 无 参考 基因 组 序列 的 样品 ,因此 
在 获得 clean reads 后 需要 进行 拼接 ,从 而 获取 后 续 
分 析 所 需 的 参考 序列 ( Grabherr et al., 2011; Hansen 
et al., 2012)。 

Trinity 利用 deBruijn 图 论 的 基本 原理 ,组 装 出 
全 长 转录 本 。 基 于 对 Trinity 组 装 出 来 的 转录 本 序 
列 , 取 每 个 基因 中 最 长 的 转录 本 序列 作为 unigene， 
































将 其 逮 住 , 放 入 扎 好 透气 孔 的 离心 管 中 ,每 个 离心 管 
塞 人 一 小 片花 椒 叶 片 。 带 回 实验 室 置 于 25Y% 下 饲 
养 , 选 取 1 000 头 活 动能 力 较 强 的 成 虫 进行 虫 体 清 
洗 之 后 在 解剖 镜 下 进行 肉 雄 鉴别 ,于 冰 上 迅速 进行 
各 组 织 分 离 .收集 ,迅速 用 液 氮 冷冻 并 保存 于 - 80Y 
低温 冰箱 中 备用 。 














再 对 unigene 长 度 分 布 进行 统计 。 将 统计 到 的 不 同 
长 度 的 unigene 分 别 与 NR，KEGG 和 COG 数据 库 
进行 序列 比 对 和 信息 分 析 ( 期 望 值 E <1.00E™™ )。 
在 此 基础 上 完成 全 部 unigene 序列 的 方向 定位 和 同 
源 序列 信息 汇总 。 将 不 同样 品 组 装 结果 中 多 个 可 变 
剪接 的 转录 本 用 Chrysalis 聚 类 得 到 unigene 库 。 通 
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过 unigene 在 Mapped Reads 上 的 分 布 , 检 验 片段 化 
的 随机 程度 、 基 因数 目的 饱和 情况 。 测 序数 据 量 可 
正 向 反映 测序 基因 的 数目 ,数据 量 越 大 得 到 的 基因 
数目 越 多 。 通 过 unigene 核酸 序列 与 数据 库 的 
BLAST 比 对 ,可 得 到 unigene 的 蛋白 功能 注释 信息 ， 
比 对 中 选择 参数 E-value 不 大 于 10 。 
1.5 ”unigene 的 功能 注释 

采用 一 些 主要 功能 数据 库 , Swiss-Prot，NR ， 
NT, PFAM, GO, KO, egeNOG, BLASTP, BLASTX 
等 对 基因 和 蛋白质 功能 注释 。 
1.6 花椒 窄 吉本 化 学 感受 相关 基因 系统 发 育 树 分 析 

利用 NCBI 在 线 软件 进行 花椒 窗 吉 丁 unigene 
的 BLAST 比 对 ,将 比 对 相似 度 较 高 的 鞘翅 目 昆虫 的 
化 学 感受 相关 基因 序列 下 载 下 来 ,利用 MEGA 软件 
构建 系统 发 育 树 , 比较 花椒 窗 吉 丁 候选 气味 受 体 
(odorant receptor，OR ) 、 离 子 型 受 体 (ionotropic 
receptor, IR) 和 味觉 受 体 ( gustatory receptor，GR ) 基 
因 与 畏 怒 目 其 他 昆虫 相关 基因 的 近 源 关 系 。 
1.7 花椒 窄 吉 丁 化 学 感受 相关 基因 的 表达 量 分 析 
通过 RPKM (reads per kilobase per millon 
mapped reads ) 值 来 衡量 基因 在 花椒 鹤 吉 丁 成 虫 不 
同 组 织 中 的 表达 水 平 ,RPKM 是 利用 RNA-Sed 技术 
用 来 定量 估计 基因 表达 值 的 有 效 工具 ,能 够 消除 基 

































































达 量 差异 。 其 计算 公式 为 : 
公式 中 民 为 唯一 比 对 到 某 个 基因 的 reads 数 ,V 
为 唯一 比 对 到 基因 的 总 reads 数 ,L 为 此 基因 的 长 度 。 
10° xR 
NILZ103 
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2 结果 


2.1 花椒 窄 吉 丁 转录 组 的 测序 和 组 装 

利用 Ilumina NovaSeq6000 平台 测序 ,组 装 得 到 
80 320 条 unigenes 和 169 398 条 contigs( 表 1),G+ 
C 比例 分 别 是 37.63% 和 39. 18% , 碱 基 总 数 分 别 是 
183 683 635 和 66 565 475 。contigs 和 unigenes 的 平 
均 长 度 分 别 为 1 084.33 和 828.75 bp( 表 2) ,N50 长 
度 分 别 为 2 108 和 1 438 bp。 长 度 在 200 ~ 400 bp 
的 unigenes 有 37 374 条 , 占 全 部 unigenes 的 
46. 53% ;而 400 ~600 bp 的 有 14 108 条 占 17.56%，; 
其 余 长 度 的 unigenes 数量 占 比 均 小 于 10% (图 1)， 
表明 测序 结果 和 组 装 效果 很 好 ,满足 进一步 开展 基 
因 功 能 分 析 要 求 。 

表 1 花椒 窄 吉 丁 转录 组 组 装 结果 统计 


Table 1 Summary of assembly for the transcriptome 








of Agrilus zanthoxylumi 








数量 (G+C)% 碱 基 总 数 
长 度 和 测序 量 差 异 对 计算 基 表 达 的 影 响 计算 Number Total number of bases 
得 到 的 基因 表达 量 可 直接 用 于 比较 不 同样 品 间 的 基 Coniss 。 169398 39.18 66 565 475 
因 表达 差异 ,此 文中 旨 在 比较 不 同 组 织 中 的 基因 表 。 一 
表 2 花椒 窄 吉 丁 转录 组 序列 长 度 统计 
Table 2 Statistics of sequence length of the transcriptome of Agrilus zanthoxylumi 
N50 N90 最 小 长 度 (bp) 最 大 长 度 (bp) 数量 平均 长 度 (bp) 
(bp) (bp) Minimum length Maximal length Number Mean length 
Contigs 2 108 399 201 28 219 169 398 1 084.33 
Unigenes 1 438 313 201 28 219 80 320 828.75 





N50: 序列 从 大 到 小 排列 , 当 其 累计 长 度 刚 刚 超过 全 部 序列 总 长 度 50% 时 ,最 后 一 个 序列 的 大 小 即 为 N50 的 大 小 The sequences are arranged from 
large to small, and when the accumulative length just exceeds 50% of the total length of all sequences, the size of the last sequence is N50; N90. 序列 
从 大 到 小 排列 , 当 其 累计 长 度 刚 刚 超过 全 部 序列 总 长 度 90% 时 ,最 后 一 个 序列 的 大 小 即 为 N90 的 大 小 The sequences are arranged from large to 


small, and when the cumulative length just exceeds 90% of the total length of all the sequences, the size of the last sequence is the size of N90. 

















Q30 值 为 99. 9% ,说 明 碱 基 识 别 率 较 高 ,测序 
质量 高 数据 可 靠 性 好 ,对 应 关系 如 表 3 所 示 。 

为 保证 分 析 准 确 性 ,以 Q30 碱 基 百 分 比 作为 指 
标 进行 统计 clean data 的 质量 ,Q30 碱 基 百分比 越 大 
说 明 测 序 错误 率 小 于 0.1% 的 碱 基 在 总 碱 基 中 的 比 
例 越 大 。 图 2 以 所 有 花椒 罕 吉 丁 不 同 组 织 ( 触 角 、 
头 \ 胸 \ 腹 、 足 过) 为 横 坐 标 ,以 Q30 百分比 作为 对 
应 的 纵 坐标 直观 地 反映 了 最 终 过 滤 后 的 数据 可 靠 























性 。 所 有 组 织 碱 基 百分比 都 在 91% 以 上 ,说 明 测 序 
数据 可 靠 可 进行 后 续 分 析 。 

2.2 花椒 窄 吉 丁 转录 组 基因 功能 注释 

通过 选择 BLAST 参数 E-value 反 le ”和 
HMMER 参数 E-value<1e-" ,最终 获得 了 80 320 个 
有 注释 信息 的 unigenes, 数据库 注释 结果 表明 ,NR 
数据 库 注释 32 470 条 (40. 43% ) ，BLASTX 注释 
29 382 条 (36.58%), GO 数据 库 注释 27 488 条 















































5 期 杨 平等 : 花椒 罕 吉 丁 转录 组 及 化 学 感受 相关 基因 的 分 析 551 











70000 


60000 


50000 


2 40000 


数量 
Number 


30000 


20000 


10000 





Fig. 1 


国 Contigs 
国 Unigenes 





序列 长 度 (bp) 
Sequence length 


图 1 花椒 窗 吉 丁 转录 组 contigs 和 unigenes 长 度 分 布 


Length distribution of contigs and unigenes in the transcriptome of Agrilus zanthoxylumi 


表 3 Inumina Casava 碱 基 识 别 与 Phred 分 值 之 间 的 简明 对 应 关系 
Table 3 The concise correspondence between lllumina Casava base recognition and Phred score 
































Phred 分 值 不 正确 的 碱 基 识 别 碱 基 正 确 识别 率 ( % ) Q 分 值 
Phred score Incorrect base-recognition Correct base recognition rate Q-score 
10 1/10 90.00 Q10 
20 1/100 99.00 Q20 
30 1/1 000 99.90 Q30 
40 1/10 000 99.99 Q40 
93.00] 表 4 花椒 窄 吉 丁 转录 组 基因 功能 注释 结果 统计 
92.501 Table 4 Statistics of unigene annotation of the 
全 92.00] transcriptome of Agrilus zanthoxylumi 
品 数据 库 Unigene 数量 百分比 
91.501 
Database Number of unigenes Percentage 
3509 Swiss-Prot 25 941 32.30 
90.50 触角 T T 胸 T 用 T 足 T 旭 NR 32 470 40.43 
y 回 | 
Antenna Head Thorax Abdomen Leg Wing Ny 12:058 Loo 
组 织 Tissues PFAM 18 038 22.46 
GO 27 488 34.22 
如 2 ”花椒 窄 吉 丁 转录 组 e 4 质量 统 j 
图 2 化 椒 罕 富 本 转录 组 clean data 的 质 统计 Re a a 
Fig. 2 Quality Sashes of clean data a the transcriptome eggNOG 18 013 22 43 
有 of Agrilus 2 请 18 921 aese 
Q30: 过 滤 后 ,总 序列 中 质量 值 大 于 30( 错误 率 小 于 0.1% ) 的 碱 BA 





基数 的 比例 Percentage of bases with a quality score greater than or 


equal to 30 (incorrect base recognition rate: <0.1% ). 


(34.22% ) ，SwissProt 库 注释 25 941 条 (32. 30% ) ， 
BLASTP 注释 18 921 条 (23. 56% ) ，PFAM 数据 库 
注释 18 038 条 (22. 46% ) ，eggNOG 注释 18 013 条 
(22. 43% ) ，KO 数据 库 注释 17 312 条 (21.55% ) 以 
及 NT 数据 库 注释 12 056 条 (15.01% )( 表 4) 。 

在 BLASTP, NR, NT 和 BLASTX 这 4 个 数据 库 
中 都 得 到 注释 的 unigene 数量 是 6 221 条 (图 3) ;在 
NR,NT 和 BLASTX 这 3 个 数据 库 注 释 的 结果 相对 


29 382 36.58 


比较 多 ,所 测 得 到 的 花椒 罕 吉 丁 的 unigene 数量 能 
在 现 阶段 各 数据 库 中 都 注释 到 一 定 的 数量 。 

通过 与 NR 数据 库 比 对 注释 ,花椒 罕 吉 丁 转录 
组 注释 到 其 他 物种 分 布 图 显示 (图 4) ,与 赤 拟 谷 盗 
Tribolium castaneum 的 相似 基因 所 占 比 例 较 高 
(26.80%)， 其 后 依次 是 小 家 上 鼠 Mus musculus 
(13. 87% ) , 性 犀 金 怨 Oryctes borbonicus (5.26% ) ， 
山 松 甲虫 Dendroctonus ponderosae (4. 10% ) ， 黑 凡 
Lasius niger(2.36% ) ,其 他 物种 占 47.61% 。 
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图 3 花椒 窗 吉 丁 转录 组 基因 功能 注释 结果 韦 恩 图 


Fig. 3 Venn diagram of the annotation results of the 





transcriptome of Agrilus zanthoxylumi 
图 中 数值 表示 注释 到 某 个 数据 库 的 花椒 窄 吉本 unigene 数量 , 交 
部 分 的 数值 则 表示 同时 注释 在 重 钱 部 分 数据 库 的 unigene 数 


量 。The numerals in the figure represent the number of unigenes 











annotated into a certain database, while the numerical value of the 
intersection part indicates the number of unigenes annotated in the 


overlapping part of the database. 


赤 拟 谷 盗 Tribolium castaneum 
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4 ”NR 数据 库 中 注释 的 花椒 罕 吉 丁 unigenes 的 
物种 分 布 


Fig. 4 Species distribution of Agrilus zanthoxylumi 

















unigenes in the NR database 


2.3 ”花椒 窄 吉 丁 unigenes 的 功能 分 类 

经 G0 数据 库 比 对 到 27 488 个 unigenes 可 分 3 
个 功能 类 别 : 分 子 功能 (molecular function ) 、 细 胞 组 
分 (cellular component ) 和 生物 学 进程 (biological 


大 分 子 相 关 的 unigenes 所 占 比例 分 别 是 82. 91% ， 
53.38% , 45.25% , 30.62% 和 30. 62% 。 而 与 分 子 
功能 相关 的 unigenes 有 14 类 ,最 多 的 是 结合 相关 
的 , 占 67. 82% , 分 子 功能 中 具有 催化 活性 的 
unigenes 上 二 48.14% 。 

将 花椒 罕 吉 本 的 unigenes 与 KOG 数据 库 进 行 
比 对 并 根据 其 功能 进行 分 类 统计 (图 6) ,注释 到 
24 614 个 unigenes 共 25 类 。 注 释 到 unigene 最 多 的 
功能 类 别 是 普通 功能 , 共 4 666 个 ,其 次 为 翻译 后 修 
饰 、 蛋 白质 折 三 和 分 子 伴侣 功能 共 2 548 个 ,而 注释 
到 最 少 的 是 与 细胞 移动 相关 , 仅 40 个 。 所 注释 的 功 
能 涉及 到 花椒 窄 吉本 的 绝 大 多 数 生 命 活动 。 

将 80 320 个 unigenes 映射 到 KEGG 数据 库 中 ， 
注释 为 代谢 通路 相关 的 序列 共有 21 104 条 , 占 总 数 
的 26.27% ,把 这 些 代谢 通路 分 为 五 大 类 ,分 别 是 细 
胞 过 程 环境 信息 处 理 过 程 、 遗 传 信息 处 理 过 程 .新 
陈 代谢 过 程 和 有 机 体系 统 过 程 ; 共 涉 及 到 35 种 代谢 
途径 ,其 中 注释 的 unigene 最 多 的 是 翻译 途径 , 共 涉 
及 到 2 037 条 unigenes, 占 18.44% ;其 次 为 信号 转 导 
途径 , 共 涉 及 到 1 910 条 unigenes, 占 17.29% ;第 3 
为 碳水 化 合 物 代谢 , 共 涉 及 到 1 336 条 unigenes, 占 
12. 09% 。 其 他 更 多 的 代谢 途径 涉及 运输 和 分 解 代 
谢 、 氨 基 酸 代谢 等 (图 7)。 

对 花椒 罕 吉 丁 转录 组 数据 查找 并 对 查找 到 的 候 
选 序列 在 NR 数据 库 中 进行 比 对 验证 后 注释 到 8 个 
气味 受 体 (odorant receptor，OR ) 基因 7 个 离子 型 受 
体 ( ionotropic receptor，IR ) 基因 和 6 个 味觉 受 
( gustatory receptor，GR ) 基因 ,共计 21 个 ( 表 5,6)。 
根据 基因 命名 规则 分 别 对 其 命名 为 AzanOR1 - 8， 
AzanlR1 -7 和 AzanGR1 -6。 

2.4 ”花椒 窄 吉 丁 转录 组 化 学 感受 相关 基因 的 系统 
发 育 分 析 

图 8 所 示 , 花椒 窄 吉 丁 的 OR 与 白蜡 窗 吉 
4sgrilus planipennis， 黄 粉 虫 Tenebrio moliior， 上 暗黑 鲁 
金龟 Holotrichia parallela， 大 猿 叶 甲 Colaphellus 
bowringi,， 大 斑 沟 死 苹 Hycleus phaleratus 和 星 天 牛 
Anoplophora chinensis 6 种 鞘翅 目 昆 虫 OR 的 系统 发 











process) 三 大 类 共 59 个 过 程 ( 图 5)。 从 图 中 结果 可 


育 树 ,从 图 可 知 AzanOR1，AzanOR2 和 AzanOR6 与 





知 参 与 生物 学 进程 的 unigenes 有 27 类 ,其 中 参与 细 
胞 进程 的 unigenes 占 总 体 的 73. 24% ,而 代谢 进程 
的 占 60. 10% ,这 两 个 过 程 所 占 比 例 最 多 ;其 次 是 生 
物 调 控 和 细胞 组 分 组 织 或 生物 合成 分 别 占 38. 22% 
和 22.71% 。 细 胞 组 分 相关 的 unigenes 有 18 类 ,其 
中 与 细胞 组 分 、 细 胞 器 、 细 胞 器 组 分 、 膜 结构 和 复杂 


白蜡 窄 吉 丁 OR 聚 类 在 一 起 ,相似 度 高 ,AzanOR3 与 
大 斑 沟 欧 蔷 OR 相似 度 比 较 高 , 而 AzanOR4 和 
AzanOR5 则 与 暗黑 鳞 金 包 OR 的 相似 度 较 高 ;图 9 
是 花椒 窗 吉 丁 IR 与 白蜡 罕 吉 丁 Agrilus planipennis， 
大 猿 叶 甲 Colaphellus bowringi， 有 暗黑 鳃 金龟 Holotrichia 
parallela, 山 松 甲虫 刀 ponderosae，, 黄粉 虫 Tenebrio 
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图 5 花椒 罕 吉 丁 unigenes CO 分 类 


Fig. 5 GO classification of Agrilus zanthoxylumi unigenes 


molitor 和 大 斑 沟 死 靖 Hycleus phaleratus IR 的 系统 发 
育 树 ,从 中 可 知 花 椒 鹤 吉 丁 IR 与 白蜡 窗 吉 丁 \ 大 斑 
沟 元 靖 、 暗 黑 鳃 金龟 和 大 猿 叶 甲 IR 的 相似 度 较 高 ; 
图 10 是 花椒 窗 吉 丁 的 GR 与 白蜡 罕 吉 丁 Agrilus 
planipennis ,大 猿 叶 甲 Colaphellus bowringi， 大 斑 沟 匹 
星 天 和 牛 Anoplophora chinensis， 
赤 拟 谷 资 Tribolium castaneum 和 铜绿 丽 金 包 








苹 Hcleus phaleratus， 





Anomala corpulentia GR 的 系统 发 育 树 ,花椒 罕 吉 

GR 与 白蜡 窄 吉 丁 、 大 猿 叶 甲 和 赤 拟 谷 盗 GR 的 相似 

度 较 高 。 这 3 类 受 体 都 与 花椒 罕 吉 丁 的 近 源 种 之 一 

的 白蜡 罕 吉 丁 GR 的 相似 度 较 高 ,在 一 定 程度 上 说 

丁 转录 组 数据 基因 的 准确 性 ,也 为 今 
这 类 基因 的 鉴定 提供 参照 和 思 
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W: Extracellular structures 





Y: Nuclear structure 
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功能 分 类 Function class 


图 6 花椒 窗 吉 丁 unigenes 的 KOG 功能 分 类 


Fig. 6 KOG function classification of Agrilus zanthoxylumi unigenes 


Z: Cytoskeleton 





内 分 泌 系 统 Circulatory system 加 天 到 日 360 
消化 系统 Digestive system IE 到 461 
老化 Aging |0 
免疫 系统 Immune systerm 300 
排泄 系统 Excretory system 辆 唆 晶 267 
循环 系统 Endocrine systerm [到 1056 
细胞 群体 Cellular community 到 到 到 到 唱 591 
神经 系统 Nervous system 天天] 593 
感觉 系统 Sensory system 国王 170 
环境 适应 性 Environmental adaptation 国 国 186 
发 育 Development 到] 296 
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其 他 氨基 酸 代谢 Metabolism of other amino acids [El396 
核 苷 酸 代谢 Nucleotide metabolism [555 
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总 览 Overview Eeea1095 
氨基 酸 代 谢 Amino acid metabolism 国 守 9 1095 
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翻译 Translation [ESaaaaaaaaaaaaaeaeal2037 
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膜 运输 Membrane transport 国 107 
信号 转 导 Signal transduction El1910 
RNA 家 族 RNA family |0 
信号 分 子 及 作用 Signaling molecules and interaction 辆 天 到 | 293 
细胞 运动 Cell motility 到 327 
细胞 生长 与 死亡 Cell growth and death 0 536 
运输 和 分 解 代 谢 Transport and catabolisrn [ER 1244 
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7 花椒 罕 吉 丁 unigenes KEGG 代谢 通路 图 
Fig. 7 KEGG metabolic pathway map of Aegrilus zanthoxylumi unigenes 
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表 5 花椒 窄 吉 丁 转录 组 化 学 感受 相关 基因 注释 结果 统计 
TableS Annotated chemoreception-related genes in the transcriptome of Agrilus zanthoxylumi 
化 学 感受 相关 基因 注释 的 基因 数量 
Chemoreception-related genes Number of annotated unigenes 
气味 受 体 基 因 Odorant receptor (OR) gene 8 
离子 型 受 体 基 因 Ionotropic receptor (IR) gene 7 
味觉 受 体 基因 Gustatory receptor (GR) gene 6 
表 6 花椒 窄 吉 丁 转录 组 化 学 感受 相关 基因 NR 中 的 注释 结果 
Table 6 Chemoreception-related genes in Agrilus zanthoxylumi transcriptome annotated in the NR database 
基因 基因 编号 NR 评分 NR 估算 值 描述 
Genes Gene ID NR score ”NR estimated value Description 
AzanOR1 TRINITY_DN29115_c0_gl 671 E:0 PREDICTED :odorant receptor coreceptor [ Tribolium castaneum | 
AzanOR2 TRINITY_DN31824_c5_gl 111 了 :4e -24 Odorant receptor OR11{ Colaphellus bowringi | 
AzanOR3 TRINITY_DN24144_c0_gl 313 了 :3e -99 Odorant receptor 14 [ Anomala corpulenta ] 
AzanOR4 TRINITY_DN30901_c0_gl 384 E:2e—118 Hypothetical protein D910_08613 [ Dendroctonus ponderosae | 
AzanORS TRINITY_DN33665_c3_gl 90.9 E:2e-16 PREDICTED : uncharacterized protein LOC103313324 [ Tribolium castaneum | 
AzanOR6 TRINITY_DN31824_c5_gl 111 了 :4e -24 Odorant receptor OR11 [ Colaphellus bowringi | 
AzanOR71 TRINITY_DN25958_cl_g3 120 了 :4e -28 Odorant receptor 14 [ Anomala corpulenta ] 
AzanOR8 TRINITY_DN25958_cl_g2 105 了 上 :8e -25 Odorant receptor 14 [ Anomala corpulenta ] 
AzanlR1 TRINITY_DN29513_cl_gl 725 E:0 Ionotropic receptor 8a[ Colaphellus bowringi | 
AzanlR2 TRINITY_DN29688_c0_gl 1 360 E:0 Chemosensory ionotropic receptor IR6 [ Tenebrio molitor | 
AzanlR3 TRINITY_DN26923_cl_g4 166 E:3e -45 Chemosensory ionotropic receptor IR3, partial [ Tenebrio molitor | 
AzanlR4 TRINITY_DN28534_c3_g2 660 E:0 Chemosensory ionotropic receptor 21a, partial [Anomala corpulenta | 
AzanlRS TRINITY_DN27218_c0_gl 325 上 :Se -96 Chemosensory ionotropic receptor IRS [ Tenebrio molitor | 
AzanlR6 TRINITY_DN28534_c3_g2 660 E:0 Chemosensory ionotropic receptor 21a, partial [Anomala corpulenta | 
AzanlR7 TRINITY_DN29688_c0_gl 1 360 E:0 Chemosensory ionotropic receptor IR6 [ Tenebrio molitor | 
4zamCR1 TRINITY_DN34348_c2_g2 66.2 E:6e -09 PREDICTED : gustatory receptor for sugar taste 43a-like [Stomoxys calcitrans | 
AzanGR2 TRINITY_DN14380_c0_gl 135 了 上 :0e -33 Gustatory receptor candidate 9[ Tribolium castaneum | 
AzanGR3 TRINITY_DN27483_c0_g3 221 E:le—64 Gustatory receptor candidate 29[ Tribolium castaneum | 
AzanGR4 TRINITY_DN29925_cl_g3 157 E:2e -39 PREDICTED : putative gustatory receptor 28b [ Tribolium castaneum | 
AzanGR5 TRINITY_DN22765_c0_gl 210 E:2e -62 Hypothetical protein YQOE_06884, partial [ Dendroctonus ponderosae] 
AzanGR6 TRINITY_DN27755_c0_g2 无 Null 无 Null 无 Null 








2.5 花椒 窄 吉 丁 转录 组 化 学 感受 相关 基因 的 表达 量 


分 析 


对 3 类 化 学 感受 


相关 基因 在 花椒 罕 吉 丁 成 虫 不 


同 组 织 中 (由 于 肉 虫 头 部 转录 组 数据 缺失 ， 
部 数据 在 此 不 加 以 分 析 和 讨论 


比 发 现 ,3 类 化 学 感受 
表达 量 相 对 其 他 组 




















著 的 为 Azan0R4 ;腹部 中 3 类 化 学 感受 
达 量 雌 、 雄 差异 最 显著 的 为 4zan0R5 ; 雌 、 雄 翅 中 3 


类 化 学 感受 


LA 














) 的 RPKM 值 进 行 
相关 基因 在 肉 、 0 
织 中 要 高 , 且 4zonrOR1 表达 量 
雌雄 差异 最 为 显著 ;胸部 中 表达 量 雌 、 雄 差异 最 显 
相关 基因 表 























相关 基因 表达 量 差 异 最 














车 的 为 


AzanIRS 和 AzanlR7 ; 而 此 、 雄 成 虫 足 中 3 类 化 学 感受 


相关 基因 表达 量 没有 差异 (图 11)。 
3 讨论 








新 一 代 高 通 量 测序 技术 的 不 断 发 展 ,已 彻底 改 
变 了 转录 组 学 的 研究 ,使 RNA-Sed 无 需 预先 设计 探 
针 即 可 对 特定 条 件 下 任何 生物 生长 发 育 阶段 整体 转 
录 活 动 进行 测序 ,准确 探测 到 各 种 条 件 下 的 基因 表 
达 情 况 ,发 现 了 许多 未 知 的 分 子 调 控 机 制 ( 刘 红 亮 
等 , 2013; 朱 志 明 等 , 2016; 王 文 龙 等 , 2017 ) ,一 
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图 8 邻接 法 构建 的 基于 氨基酸 序列 的 花椒 窄 吉 丁 和 其 他 物种 OR 系统 发 育 树 





Fig. 8 Phylogenetic tree of ORs from Agrilus zanthoxylumi and other insect species by 


using neighbor-joining method based on the amino acid sequence 


蛋白 来 源 物种 Origin species of proteins: AplanOR: 白蜡 窄 吉 丁 Agrilus planipennis; TmolOR: 黄粉 虫 Tenebrio molitor; HparOR: 暗黑 鲁 金 
Colaphellus bowringi; HphaOR: 大 斑 沟 
化 树 中 序列 名 为 GenBank 登录 号 + 蛋白 名 。 以 bootstrap =1 000 构 建 进 














Holotrichia parallela; CbowOR: 大 猿 叶 
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Hycleus phaleratus; AchiOR: 星 天 牛 Anoplophora chinensis. 进 














化 树 。The names of sequences in the phylogenetic tree are GenBank 


accession no. plus protein name. Bootstrap support values presented as percentages are based on 1 000 replicates. 
Pp P: Pp supp Pp Pp 8 Pp 


重要 功能 基因 已 被 挖掘 并 加 以 利用 (Sangwan et al.， 
2013 ; Campos-Parra et al., 2016; Xu et al., 2016; 
Nguyen et al., 2017 ) 。 花 椒 窗 吉 丁 转录 组 测序 、 组 
装 和 分 析 后 共 得 到 80 320 条 基因 ,169 398 条 转录 
本 ;contigs 和 unigenes 的 G +C 比例 分 别 是 39. 18% 
和 37. 63%; 所 有 样本 的 Q30 碱 基 比 例 不 低 于 
91. 00% ,unigenes 的 N50 长 度 为 1438 bp; 为 保证 测 
序 质量 Q30 碱 基 比 例 不 得 低 于 80% ,而 N50 的 值 越 
大 表示 序列 拼接 的 完整 性 越 好 ( 王 兴 春 等 , 2015 ) ， 
说 明 本 实验 获得 的 测序 数据 质量 好 ,保证 了 转录 组 
分 析 的 准确 性 ,其 中 测序 错误 率 与 碱 基质 量 有 关 , 受 
测序 仪 本 身 、 测 序 试剂 .样品 等 多 个 因素 共同 影响 。 
在 NR 数据 库 中 花椒 窄 吉 丁 转录 组 注释 到 其 他 的 物 


























种 中 赤 拟 谷 盗 所 占 比例 为 26. 80% ,其 后 依次 是 小 
家 鼠 13.87% , 星 犀 金 包 5.26% , 山 松 甲虫 4.10%， 
黑 蚁 2.36% ,其 他 物种 占 比 例 最 高 (47. 61% )。 在 
G0 数据 库 中 的 比 对 到 27 488 个 unigenes, 可 分 为 分 
子 功能 、 细 胞 组 分 和 生物 过 程 三 大 类 共 59 个 过 程 ， 
但 在 G0 数据 库 中 注释 到 unigenes 仅 占 总 unigenes 
数 的 34.22% ,未 得 到 注释 的 这 些 可 能 是 因为 新 的 
序列 ,未 在 数据 库 中 存在 或 者 是 进行 非 编 码 的 ;KOC 
数据 库 功 能 注释 的 unigene 根据 其 功能 分 为 25 类 ， 
注释 到 unigene 最 多 功能 类 别 的 是 普通 功能 , 共 
4 666 个 ,其 次 为 翻译 后 修饰 蛋白 质 折 释 和 分 子 伴侣 
功能 共 2 548 个 ,而 注释 到 最 少 的 是 与 细胞 移动 相 
关 , 仅 40 个 ;将 这 80 320 个 unigenes 映射 到 KEGG 
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Fig. 9 Phylogenetic tree of IRs from Agrilus zanthoxylumi and other insect species by 


using neighbor-joining method based on the amino acid sequence 
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IR 来 源 物种 Source species of IRs: AplanIR: 白蜡 罕 吉 丁 Agrilus plan 





Holotrichia parallela; DponIR: 山 松 甲 上 
树 中 序列 名 为 GenBank 登录 号 + 蛋白 名 。 以 bootstrap =1 000 构 如 


no. plus protein name. Bootstrap support values presented as percentages a 
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数据 库 中 ,注释 到 代谢 通路 35 个 ,其 相关 的 序列 共 
有 21 104 条 , 占 总 数 的 26.27% ,把 这 些 代谢 通路 分 
为 五 大 类 ,分别 是 细胞 过 程 、 环 境 信息 处 理 过 程 . 遗 
传 信息 处 理 过 程 新 陈 代谢 过 程 和 有 机 体系 统 过 程 ， 
这 些 代谢 通路 的 分 类 与 已 报道 的 研究 中 并 没什么 不 
同 。 本 研究 对 花椒 窄 吉 丁 转录 组 原始 数据 的 公开 及 
分 析 ,为 种 属 鉴定 及 遗传 多 样 性 的 分 析 提 供 了 基础 
资料 ,同时 也 为 了 解 其 基因 功能 .生物 学 特征 以 及 相 
关 的 生理 学 研究 提供 数据 参考 。 

本 研究 首次 从 花椒 罕 吉 丁 转录 组 数据 中 挖掘 到 
8 个 气味 受 体 基因 7 个 离子 型 受 体 基因 和 6 个 味觉 
受 体 基因 ,共计 21 个 。 味 觉 受 体 主要 包含 二 氧化 碳 


























R Derdroclomas pomderosae; TmolIR : 黄粉 
进化 树 。The names of sequences in the phylogenetic tree are GenBank accession 
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ipennis; CbowIR: 大 猿 叶 甲 Colaphellus bowringi; HparIR: 暗黑 鳗 金 
R Tenebrzo molitor; HphaIR: 大 斑 沟 欧普 Hycleus phaleratus. 进 




















re based on 1 000 replicates. 


检测 受 体 (Kwon et al., 2006) ,然而 ,大 多 数 GRs 在 
味觉 器 官 的 味觉 受 体 神经 元 中 表达 并 参与 接触 化 学 
感受 ,这 些 GRs 通常 检测 不 同 的 糖 、 藻 味 化 合 物 和 
接触 信息 素 (Robertson and Kent，2009 ) 。 昆 虫 ORs 
是 含 7 个 跨 膜 结构 域 的 蛋白 (Clyne ei al., 1999) ,与 
GRs 是 同一 个 受 体 基因 超 家 族 的 (Vosshall et al.， 
1999 ) ,但 ORs 和 GRs 相互 之 间 的 序列 相关 性 较 小 
( Robertson et al., 2003 ) 。 而 IRs 是 一 类 参与 化 学 感 
受 的 典型 受 体 (Benton et ol.，2009) ,它们 主要 在 突 


I 
触 通讯 中 起 作用 。 不 同 昆虫 具有 不 同 数目 的 气味 接 
受 相关 基因 , 与 花椒 罕 吉 丁 近 源 的 苹果 小 吉 丁 4. 
mali 中 化 学 感受 相关 基因 分 别 有 17 条 OR 基因 、17 
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图 10 ”邻接 法 构建 的 基于 氨基 酸 序列 的 花椒 罕 吉 丁 和 其 他 物种 GR 系统 发 育 树 
Fig. 10 Phylogenetic tree of GRs from Agrilus zanthoxylumi and other insect species by 


using neighbor-joining method based on the amino acid sequence 


GR 来 源 物种 Source species o 
Hycleus phaleratus; AchiGR: 星 天 4 
树 中 序列 名 为 GenBank 登录 号 + 蛋白 名 。 以 bootstrap =1 000 构建 进 























GRs: AplanGR: 白蜡 窄 吉 丁 Agrilus planipennis; CbowGR: 大 猿 叶 甲 Colaphellus bowringi; HphaGR: 大 斑 沟 芜 靖 
F Anoplophora chinensis; TeasGR : 赤 拟 谷 盗 Tribolium castaneum; AcorGR; 铜绿 丽 金 包 4Anomala corpulenta. 进化 


EE 进化 树 。The names of sequences in the phylogenetic tree are GenBank accession 








no. plus protein name. Bootstrap support values presented as percentages are based on replicates. 
plus prot Bootstrap support values p ted as p' tag' based on 1 000 replicat 


条 IR 基因 和 9 条 GR 基因 ( 崔 晓 宁 , 2018) ,但 其 中 
花椒 罕 吉 丁 OR 和 IR 基因 的 数量 与 苹果 小 吉 丁 的 
数量 相差 较 大 ,可 能 是 因为 在 花椒 罕 吉 仅 危 害 花椒 
在 接受 气味 分 子 上 存在 一 定 特异 性 ,化 学 感受 相关 
基因 的 数量 上 可 能 会 有 所 减少 。 因 此 有 关 花 椒 罕 吉 
丁 化 学 感受 相关 基因 的 组 织 表 达 特 异性 、 功 能 和 接 


受气 味 分 子 的 特异 性 还 需 进一步 实验 分 析 验 证 。 
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Chemoreception-related genes 


图 11 花椒 罕 吉 丁 化 学 感受 相关 基因 在 成 虫 不 同 组 织 中 的 RPKM 值 
Fig. 11 The RPKM values of chemoreception-related genes in different tissues of Agrilus zanthoxylumi adults 
R: 气味 受 体 Odorant receptor; IR: 离子 型 受 体 Ionotropic receptor; GR: 味觉 受 体 Gustatory receptor. A: 触角 Antenna; B: 胸 Thorax; C: 腹 
Abdomen; D: 足 Leg; E: 翅 Wing. 本 研究 采用 RPKM (每 百 万 reads 中 来 自 某 一 基因 每 千 碱 基 长 度 的 reads 数 ) 来 表示 基因 的 表达 量 , 即 图 中 
纵 坐 标 RPKM 值 表示 基因 的 相对 表达 量 。In this study, the number of reads per million from a certain gene per 1 000 base length ( RPKM) was used 


to represent the gene expression level, that is, the ordinate RPKM value in the figure represents the relative expression level of the gene. 
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